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In neutral cold quark matter that is sufficiently dense that the strange quark mass Ms is unim- 
portant, all nine quarks (three colors; three flavors) pair in a color-flavor locked (CFL) pattern, and 
all fermionic quasiparticles have a gap. We recently argued that the next phase down in density 
(as a function of decreasing quark chemical potential ^ or increasing strange quark mass Ms) is the 
new "gapless CFL" ("gCFL") phase in which only seven quasiparticles have a gap, while there are 
gapless quasiparticles described by two dispersion relations at three momenta. There is a continuous 
quantum phase transition from CFL to gCFL quark matter at Mf //i ~ 2A, with A the gap param- 
eter. Gapless CFL, like CFL, leaves unbroken a linear combination Q of electric and color charges, 
but it is a Q-conductor with gapless Q-charged quasiparticles and a nonzero electron density. In 
this paper, we evaluate the gapless CFL phase, in several senses. We present the details underlying 
our earlier work which showed how this phase arises. We display all nine quasiparticle dispersion 
relations in full detail. Using a general pairing ansatz that only neglects effects that are known to 
be small, we perform a comparison of the free energies of the gCFL, CFL, 2SC, gapless 2SC, and 
2SCus phases. We conclude that as density drops, making the CFL phase less favored, the gCFL 
phase is the next spatially uniform quark matter phase to occur. A mixed phase made of colored 
components would have lower free energy if color were a global symmetry, but in QCD such a mixed 
phase is penalized severely. 

PACS numbers: 



I. INTRODUCTION 

Because QCD is asymptotically free, we expect that 
matter at sufficiently high densities and/or temperatures 
will consist of almost-free quarks and gluons. However, 
over the last few years it has become clear that there is a 
rich and varied landscape of phases lying between these 
asymptotic regimes and the familiar hadronic phase at 
low temperature and density. In the region where the 
temperature is low and the density is high enough that 
hadrons are crushed into quark matter, there is a whole 
family of "color superconducting" phases Q . The essence 
of color superconductivity is quark pairing, driven by the 
BCS mechanism, which operates when there exists an 
attractive interaction between fermions at a Fermi sur- 
face. The QCD quark-quark interaction is strong, and 
is attractive in many channels, so we expect cold dense 
quark matter to generically exhibit color superconductiv- 
ity. Moreover, quarks, unlike electrons, have color and 
flavor as well as spin degrees of freedom, so many dif- 
ferent patterns of pairing are possible. This leads us to 
expect a rich phase structure in matter beyond nuclear 
density. 

Color superconducting quark matter may well occur 
naturally in the universe, in the cold dense cores of com- 
pact ("neutron") stars, where densities are above nuclear 
density, and temperatures are of the order of tens of keV. 
In future low-energy heavy ion colliders, such as the Com- 
pressed Baryonic Matter Experiment at the future accel- 
erator facility at GSI Darmstadt |3] , it could conceivably 
be possible to create color superconducting quark matter 



(or perhaps hot dense matter that is in the quark-gluon 
plasma phase but which exhibits fluctuations that are 
precursors of color superconductivity [3j). 

It is by now well-established that at asymptotic den- 
sities, where the up, down and strange quarks can be 
treated on an equal footing and the potentially disruptive 
strange quark mass can can be neglected, quark matter is 
in the color-flavor locked (CFL) phase, in which quarks of 
all three colors and all three flavors form Cooper pairs |j] . 
However, just as RHIC is teaching us about the proper- 
ties of the hot but far from asymptotically hot quark- 
gluon plasma |^, we should expect that if neutron star 
cores are made of color superconducting quark matter, 
they may not reach the densities at which CFL predom- 
inates. In this paper, as in Ref. [g, we ask what form of 
color superconducting quark matter is the "next phase 
down in density" . That is, we imagine beginning in the 
CFL phase at asymptotic density, reducing the density, 
and assume that CFL pairing is disrupted by the heav- 
iness of the strange quark before color superconducting 
quark matter is superseded by the hadronic phase. Upon 
making this assumption, we ask what form the disrup- 
tion takes and what are the properties of the resulting 
phase of dense, but not asymptotically dense, matter. 

To describe quark matter as may exist in the cores of 
compact stars, we consider quark chemical potentials /i 
of order 500 MeV at most. The strange quark mass Ms 
must then be included: it is expected to be density de- 
pendent, lying between the current mass ~ 100 MeV and 
the vacuum constituent quark mass '^ 500 MeV. In bulk 
matter, as is relevant for compact stars where we are 



interested in kilometer-scale volumes, we must further- 
more require electromagnetic and color neutrality [3, Q 
(possibly via mixing of oppositely-charged phases) and 
allow for equilibration under the weak interaction. All 
these factors work to pull apart the Fermi momenta of the 
different quark species, imposing an energy cost on the 
cross-species pairing that characterizes color-flavor lock- 
ing. At the highest densities we expect CFL pairing, but 
as the density decreases the combination of nonzero Ms 
and the constraints of neutrality put greater and greater 
stress on cross-species pairing, and we expect transitions 
to other pairing patterns. 

In this paper we study the first of these transitions, 
and work exclusively at zero temperature, which is a rea- 
sonable approximation in the interior of a neutron star 
that is more than a few seconds old. (Nonzero tempera- 
ture adds interesting new facets to the analysis 9J, that 
we shall further analyze elsewhere.) We argue that the 
CFL phase will first give way, via a continuous phase 
transition, to a new phase with gapless fermions that 
we call gapless CFL (gCFL). The transition occurs when 
Ms /fJ, — 2 A, where A is the pairing gap parameter. The 
gCFL phase has gapless modes and nonzero electron den- 
sity. Although it has the same symmetries as the CFL 
phase, gapless CFL matter is a conductor whereas CFL 
quark matter is a dielectric insulator. 



A. Summary of the CFL phase 

To set the stage for our analysis^we briefly summarize 
the properties of the CFL phase |j]. If we set all three 
quark masses to zero, the diquark condensate in the CFL 
phase spontaneously breaks the full symmetry group of 
QCD, 

[5(7(3)color]x 5[/(3)l X SUi3)HXU{l)B 



^ [c^(i)q] 

SU{3)c+l+rxZ2 



(1) 



where SU(3)co\ot and electromagnetism U{1)q are 
gauged, and the unbroken SU{3)c+l+r subgroup con- 
sists of flavor rotations of the left and right quarks with 
equal and opposite color rotations, and contains an un- 
broken gauged "rotated electromagnetism" C^(1)q [3.riO|. 
The CFL phase has the largest possible unbroken sym- 
metry consistent with diquark condensation, achieved by 
having all nine quarks participate equally in the pairing, 
and this gives the maximal pairing free energy benefit. 
Not surprisingly, ab initio calculations valid at asymp- 
totic densities confirm that the CFL phase is the ground 
state of QCD in the high density limit |ll[lll|. 

In the limit of three massless quarks described above 
there are 17 broken symmetry generators in the CFL 
phase, 8 of which become longitudinal components of 



massive gauge bosons and 9 of which remain as Gold- 
stone bosons. However, in the real world there are two 
light quark flavors, the up (u) and down (d), with masses 
< 10 MeV, and a medium- weight flavor, the strange 
(s) quark, with mass > 100 MeV. The strange quark 
therefore plays a crucial role in the phases of QCD. 
In the presence of quark masses, the eight Goldstone 
bosons coming from the breaking of chiral symmetries 
acquire masses |l|, |j, [T3| , and furthermore the CFL con- 
densate may rotate within the manifold describing these 
mesons |13j . In analyzing the response of the CFL phase 
to the strange quark mass, we shall be concerned with 
the dispersion relations describing its fermionic quasipar- 
ticles, as they signal an instability corresponding to the 
disruption of pairing itself. In this analysis, we shall ne- 
glect flavor rotations of the CFL condensate, as the di- 
rect effects of such meson condensates on the stability or 
instability with respect to pair breaking is minimal. (Me- 
son condensates would play an important indirect role if 
they were charged, but the favored meson condensation 
channels are neutral assuming that the neutrino density 
is negligible ^^■) The ninth Goldstone boson, that cor- 
responding to the spontaneous breaking of U{1)b and 
hence to superfluidity, remains massless even once quark 
masses are taken into account and therefore plays a cru- 
cial role in many low energy properties of the CFL phase, 
for example in its viscosity ll4l|, specific heat, neutrino 
opacity, and neutrino emissivity at low temperatures jl^ . 



B. (Gapless) CFL Pairing Ansatz 

To study the response of the CFL phase to a non- 
negligible strange quark mass, we use the pairing 
ansatz |j| 

(CCTsV'f) - A,e''''habi+A2e''''^eab2+A3e"^\at3 (2) 

Here i/)" is a quark of color a = (r, g, b) and flavor 
a = {u, d, s); the condensate is a Lorentz scalar, antisym- 
metric in Dirac indices, antisymmetric in color (the chan- 
nel with the strongest attraction between quarks), and 
consequently antisymmetric in flavor. The gap parame- 
ters Ai, A2 and A3 describe down-strange, up-strange 
and up-down Cooper pairs, respectively. They describe 
a 9 X 9 matrix in color- flavor space that, in the basis 
(ru, gd, bs, rd, gu, rs, bu, gs, bd), takes the form 
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We see that {rd,gu), {bu,rs) and (gs,bd) quarks pair 
with gap parameters Ai, A2 and A3 respectively, while 
the {ru, gd, bs) quarks pair among each other involving 
all the A's. The most important physics that we are 
leaving out by making this ansatz is pairing in which 
the Cooper pairs are symmetric in color, and therefore 
also in flavor. Diquark condensates of this form break no 
new symmetries, and therefore must arise in the CFL 
phase jjp J^]. However because the QCD interaction 
is repulsive between quarks that are symmetric in color 
these condensates are numerically insignificant H, 0, [13 ■ 
To find which phases occur in realistic quark matter, we 
must take into account the strange quark mass and equi- 
libration under the weak interaction, and impose neu- 
trality under the color and electromagnetic gauge sym- 
metries. The arguments that favor Q are unaffected 
by these considerations, but there is no reason for the 
gap parameters to be equal once Mg ^ 0. Much previ- 
ous work H El IHIll El 113 compared color-flavor- 
locked (CFL) phase (favored in the limit Mg ^ or 
/i -^ 00), the two-flavor (2SC) phase (favored in the limit 
Ms -^ 00), and unpaired quark matter. We gave a model- 
independent argument in Ref. [g, however, that when 
the CFL phase is disrupted, it cannot give way to either 
2SC or unpaired quark matter. Above a critical M^ / fi, 
we found that the CFL phase is replaced by a new gap- 
less CFL (gCFL) phase, not by 2SC quark matter. The 
defining (and eponymous) properties of the gapless CFL 
phase arise in its dispersion relations, not in its pattern 
of gap parameters. However, it is useful for orientation 
to list the patterns of gap parameters for all the phases 
we shall discuss: 



A3~A2 = Ai = Aci^L 


CFL 


(4) 


A3>0, Ai = A2 = 


(gapless) 2SC 


(5) 


A2 > 0, Ai = A3 = 


2SCus 


(6) 


A3 > A2 > Ai > 


gapless CFL . 


(7) 



The 2SCus phase, which was introduced in Ref. S\, must 
be analyzed for completeness because it and the 2SC 
phase have the same free energy when Mg — 0, and to 
leading order in Mg if their respective nonzero gap pa- 
rameters have the same value |3|- However, we shall show 
in Section IHI El that the 2SCus phase is never favored, 
and never gapless. 

In the remainder of this paper we construct the free 
energies and solve the gap equations for the CFL, gap- 
less CFL, 2SC, gapless 2SC 21, 22], and 2SCus phases in 
an NJL model. We show in detail how the CFL— >gCFL 
transition occurs and detail the properties of the gCFL 
phase. The gCFL phase is a Q-conductor with a nonzero 
electron density, and these electrons and the gapless 
quark quasiparticles make the low energy effective the- 
ory of the gapless CFL phase and, consequently, its as- 
trophysical properties qualitatively different from that of 
the CFL phase, even though its U{1) symmetries are 
the same. Both gapless quasiparticles have quadratic 
dispersion relations at the quantum critical point. For 



values of Af^//i above the quantum critical point, one 
branch has conventional linear dispersion relations while 
the other branch remains quadratic, up to tiny correc- 
tions. In order to evaluate the range of Mg/fi above 
the critical point within which the gCFL phase remains 
favored, we construct the 2SC and 2SCus phases and re- 
produce the 2SC^g2SC transition of Refs.[2ll|22, here 
in neutral 3-ffavor quark matter, and show that in this 
context gCFL has a lower free energy than (g)2SC(us). 
We do not complete the study of mixed phase alterna- 
tives, but we do eliminate all the most straightforward 
possibilities everywhere in the gCFL regime in M"^ / jjl ex- 
cept very close to its upper end, where gCFL, g2SC and 
unpaired quark matter have comparable free energies. At 
such large values of M^/^, however, our pairing ansatz 
is not sufficiently general to describe all the possibilities, 
as we discuss in the concluding section of this paper. Be- 
fore turning to the model analysis, which we detail in 
Section ^1 a-^d whose results we present in Section IIIII 
we conclude this introduction with a model-independent 
discussion of color and electric neutrality in QCD and 
with the model-independent argument of Ref. 



C. Color and Electric Neutrality in QCD 

Stable bulk matter must be neutral under all gauged 
charges, whether they are spontaneously broken or not. 
Otherwise, the net charge density would create large elec- 
tric fields, making the energy non-extensive. In the case 
of the electromagnetic gauge symmetry, this simply re- 
quires zero charge density. In the case of the color gauge 
symmetry, the formal requirement is that a chunk of 
quark matter should be a color singlet, i.e., its wave- 
function should be invariant under a general color gauge 
transformation. Color neutrality, meaning equality in the 
numbers of red, green, and blue quarks, is a less strin- 
gent constraint. A color singlet state is also color neutral, 
whereas the opposite is not necessarily true. However it 
has been shown that the projection of a color neutral 
state onto a color singlet state costs no extra free energy 
in the thermodynamic limit |23l |. Analyzing the conse- 
quences of the requirement of color neutrality therefore 
suffices for our purposes. 

In nature, electric and color neutrality are enforced 
by the dynamics of the electromagnetic and QCD gauge 
fields, whose zeroth components serve as chemical poten- 
tials which take on values that enforce neutrality ^ 2iM ■ 
Since we are limiting ourselves to color neutrality and not 
color singletness we have to consider only the U{l)xU{l) 
diagonal subgroup of the color gauge group. This sub- 
group is generated by the diagonal generators T3 = 
diag(i,-i,0) and Tg = diag(i,i,-|) of the SU{3) 
gauge group. Electromagnetism is generated by Q = 
diag(|, — -J, — -j) in flavor space (u, d, s). The zeroth com- 
ponents of the respective gauge fields serve as chemical 
potentials ^3 and fis coupled to T3 and Tg charges, and as 
an electrostatic potential /ig coupled to the negative elec- 



trie charge Q. (We make this last ehoice so that /ie > 
corresponds to a density of electrons, not positrons.) The 
dynamics of the gauge potentials then require that the 
charge densities, which are the derivatives of the free en- 
ergy with respect to the chemical potentials, must vanish: 



Q = 

T8 = 



dVl 

dUe 

on 

dVL 



(8) 



A generic diquark condensate will be neither electri- 
cally nor color neutral, so it will spontaneously break 
these gauge symmetries. However it may be neutral un- 
der a linear combination of electromagnetism and color. 
Indeed, any condensate of the form |2Jl is neutral with 
respect to the "rotated electromagnetism" generated by 
Q = Q — T3— iTg, so U{1)q is never broken. This means 

that the corresponding gauge boson (the "Q-photon"), a 
mixture of the ordinary photon and one of the gluons, 
remains massless. In both the CFL and gCFL phases, 
the rest of the S'C/(3)coior x U{1)q gauge group is spon- 
taneously broken, meaning that the combination of the 
photon and gluons orthogonal to the Q-photon, and all 
the other gluons, become massive by the Higgs mecha- 
nism. 

In an NJL model with fermions but no gauge fields, as 
we shall employ after pursuing model-independent argu- 
ments as far as we can, one has to introduce the chemical 
potentials /Xe, IJ.3 and ns "by hand" in order to enforce 
color and electric neutrality in the same way that gauge 
field dynamics does in QCD Q. 



D. Where does CFL pairing become unstable? 

We conclude this introduction with the model- 
independent argument of Ref. 6] that determines the 
density at which the CFL phase becomes unstable. The 
gap equations for the three A's will turn out to be cou- 
pled, but we can, for example, analyze the effect of a 
specified Ai on the gs and bd quarks without reference 
to the other quarks. It turns out that gs-bd pairing is the 
first to break down, and this instability is what catalyzes 
the CFL^gCFL transition. 

The leading effect of Mg is like a shift in the chemical 
potential of the strange quarks, so the bd and gs quarks 
feel "effective chemical potentials" /i^^ = 11 ~ |/i8 and 

l^fs = M + 5M8 — -2^- In the CFL phase, color neu- 
trality requires /^g = —M^/2^, a result that is model- 
independent to leading order in Mf/fj,'^ |3, [13 ■ This 
result can be understood as arising because CFL pair- 
ing itself enforces equality in the number of rd and gu 
quarks, in the number of bu and rs quarks, and in the 
number of gs and bd quarks |25l | , but in order to achieve 



neutrality the number density of (rd, gu) quarks must 
be reduced relative to that of the (bu, rs) and (gs, bd) 
quarks, and this requires a negative /ig. Because of the 
negative /xg, /xgf — ^^f — M^/fi in the CFL phase. The 
CFL phase will be stable as long as the pairing makes it 
energetically favorable to maintain equality of the bd and 
gs Fermi momenta, despite their differing effective chem- 
ical potentials f2S| . It becomes unstable when the energy 
gained from turning a gs quark near the common Fermi 
momentum into a bd quark (namely M"^ / jjl) exceeds the 
cost in lost pairing energy 2Ai. Hence, the CFL phase 
is stable when [g 



Ml 



< 2AcFL 



(9) 



For lower density, i.e. larger Af^//i, the CFL phase must 
be replaced by some new phase with unpaired bd quarks. 
One might naively expect this phase to be either neutral 
unpaired quark matter or neutral 2SC quark matter, but 
it is known that these have higher free energy than CFL 
for M'^/fj, < 4AcFL |M H^l' so this new phase, which 
must have the same free energy as CFL at the critical 
M'^/fi = 2AcFL, must be something else. In view of its 
properties that are discussed in detail in Section Hill we 
caU it gapless CFL (gCFL). 



II. MODEL AND APPROXIMATIONS 

We are interested in physics at non-asymptotic den- 
sities, and therefore cannot use weak-coupling methods. 
We are interested in physics at zero temperature and high 
density, at which the fermion sign problem is acute and 
the current methods of lattice QCD can therefore not 
be employed. For this reason, we need to introduce a 
model in which the interaction between quarks is simpli- 
fied, while still respecting the symmetries of QCD, and 
in which the effects of Mj, /ig, Ma and /xg on CFL pairing 
can all be investigated. The natural choice is to model 
the interactions between quarks using a point-like four- 
fermion interaction, which we shall take to have the quan- 
tum numbers of single-gluon exchange. We work in Eu- 
clidean space. Our partition function Z and free energy 
density $7 are then defined by 



Z = e-^^^^ ^Af fv^pVi/jexpff 



C{x) d*x 



C{x) = V(i^- 



M)^--G(v;r;^V)(^r^V') 



(10) 
where the fields live in a box of volume V and Euclidean 
time length (3 = 1/T, and ^ — fJ^j^- The interaction ver- 
tex has the color, flavor, and spin structure of the QCD 
giuon-quark coupling, T^ — 7^T . The mass matrix 
M = diag(0, 0, Mg) in flavor space. The chemical po- 
tential /I is a diagonal color-flavor matrix depending on 
/i, fj,e, fJ.3 and fis- The normalization of the four- fermion 



coupling 3G/S is as in the first paper in Ref. 1]. In real 
QCD the ultraviolet modes decouple because of asymp- 
totic freedom, but in the NJL model we have to add this 
feature by hand, through a UV momentum cutoff A in 
the momentum integrals. The model therefore has two 
parameters, the four-fermion coupling G and the three- 
momentum cutoff A, but it is more useful to parameterize 
the interaction in terms of a physical quantity, namely 
the CFL gap parameter at Ms = at a reference chemi- 
cal potential that we shall take to be 500 MeV. We shall 
call this reference gap Ag. We have checked that if we 
vary the cutoff A by 20% while simultaneously varying 
the bare coupling G so as to keep Aq fixed, then our re- 
sults change by a few percent at most. All the results 
that we present are for A = 800 MeV and for a coupling 
strength chosen such that Aq — 25 MeV. 

We now sketch the derivation of the free energy H. ob- 
tained from the Lagrangian IjlOl) upon making the ansatz 
(0) for the diquark condensate and working in the mean 
field approximation. More sophisticated derivations ex- 
ist in the literature [l|, but since we are assuming that 
the only condensate is of the form Q we simply Fierz 
transform the interaction to yield products of terms that 
appear in (j^J, and discard all the other terms that arise 
in the Fierz transformed interaction which would anyway 
vanish after making the mean field approximation. This 
yields 



G 



'^n.t = -Y.{i^PnV){i^^Pn^). 



(11) 



where 



{P,)t-c^^^ 



a[3rj 



ciji; 



(no sum over 77) (12) 



and Pri = 74^i74. The index rj labels the pairing channel: 
r] = 1,2, and 3 correspond to d-s pairing, u-s pairing, and 
u-d pairing. The overall coefficient in Hll|l is the product 
of the 3G/8 in ^ and factors of-1, 4/3, and -1/2 from 
Fierz transformations in Dirac, color and flavor space, 
respectively. 

Next, for each channel we introduce a complex scalar 
field (j)ri whose expectation value will be A,,, the strength 
of the pairing in the 77 channel, and bosonize the four- 
fermion interaction via a Hubbard-Stratonovich transfor- 
mation. The interaction Lagrangian then becomes 



Ant = ^{i>Pr,1p'^)(l)r, 



i0;(^^p,^)-:^ 



(13) 



where here and henceforth repeated tj's are summed and 
where it is understood that we are now integrating over 
the (firi as well as tp and ijj in the functional integral 
(|10(l . The functional integral is now quadratic in the 
quark fields, so the fermionic function integral can be per- 
formed. Since there are terms in the action that can vio- 
late quark number, we must use Nambu-Gorkov spinors 



* = 



_tp{p) 



* = 



(ijip) V^(-p)) (14) 



and the full Lagrange density becomes 
C 



1-5-1 
2 T 



,^v-n 



G 



where the inverse full propagator is 



s-\,)4*^i:.'' 



P, 



riVq 



P 



v'r-q 



{i>-^ + Mf 



(15) 



(16) 



We now integrate over the fermionic fields to obtain the 
effective potential for the scalar fields. We also make the 
mean-field approximation, neglecting fluctuations in the 
scalar flelds and setting 0,, to its expectation value A,,. 
The result is 



Det 



S ^{iujn,p) 



T 



1/2 



exp 



T G 



(17) 



and hence 



A„A^ 
G 



(18) 
where ujn = (2n — 1)t:T are the Matsubara frequencies. 
We do the Matsubara summation using the identity 

TY.ln{^^^)=\e\+2THl + e-\^\'^) . (19) 



In the limit of zero temperature only the first term from 
the right hand side survives, leading to the result 



n = - 



^r^^si 



47r2 
1 



iip)\dp 



(20) 



+ -(A^ + A^ + A^) 



127r2 



where the electron contribution is included, and ej{p) are 
the dispersion relations of the quasiquarks, i.e. the values 
of the energy at which the propagator diverges: 



detS''^{iej{p),p) = . 



(21) 



S~^ is a 72 X 72 matrix, but because what occurs in 
the identity (|19|1 is the combination w^ + £^, the sum 
in H20|) is understood to run over 36 roots. (This can 
be seen as removing the doubling of degrees of freedom 
introduced by using the Nambu-Gorkov formalism.) In 
the specific cases where our general ansatz becomes 2SC 
or CFL pairing, our expression (|20|l for the free energy, 
and in particular the coefficient of the A^ term, agrees 
with the expressions obtained by other methods [ij that 
do not involve Fierz transformations. 

In our numerical evaluation, we omit the antiparticle 
modes: exciting them costs of order 2/^ and they there- 
fore do not play an important role in the physics. This is 
discussed in more detail below. Neglecting the antipar- 
ticles leaves us with only 18 roots of H21|l to sum over in 
(|20|l . These correspond to 9 different dispersion relations 



describing the quasiparticlcs of differing color and flavor, 
each doubly degenerate due to spin. 

A stable, neutral phase must minimize the free energy 
(I20|l with respect to variation of the three gap parameters 
Ai, A2, A3, meaning it must satisfy 






= 0, 






= 0, 






= 



(22) 



and it must satisfy the three neutrality conditions ^. 
The gap equations (|22|) and neutrality equations © form 
a system of six coupled integral equations with unknowns 
the three gap parameters and /is, /is and fig- 

We must now find the dispersion relations £j{p), de- 
termined by the zeroes of detS'"^ which is specified by 
(|21|l . 116(1 with (pr, replaced by A,,, and (|12|l . then eval- 
uate the free energy H. using (|20ll . and then solve the 
six simultaneous equations ((HJ and ((22|l . Before carrying 
this calculation through, however, we first make a num- 
ber of simplifying approximations within the expression 
for det5-i. 

1. We neglect contributions to the condensate that are 
symmetric in color and flavor: these are known to 
be present and small 0, |3, [la| • 

2. We treat the up and down quarks as mass- 
less, which is a legitimate approximation in the 
high density regime, and we treat the constituent 
strange quark mass Mg as a parameter, rather than 
solving for an (ss) condensate. The latter approx- 
imation should be improved upon, along the lines 
ofRef. [13. 

3. We incorporate Ms only via its leading effect, 
namely as a shift — Af^/2/x in the chemical poten- 
tial for the strange quarks. This approximation ne- 
glects the difference between the strange and light 
quark Fermi velocities, whose effects arc known in 
other contexts to be small ^^. The approxima- 
tion is controlled by the smallness of M"^ / jj? . For 
this reason, in all the results that we plot we shall 
work at /i = 500 MeV and choose a coupling such 
that the CFL gap at M^ = is Aq = 25 MeV. 
We expect the CFL pairing to break down near 
M^ « 2/xAo, and choosing Aq — /i/20 ensures that 
this occurs where M'^/iJ? ~ 1/10, meaning that 
we can trust our results well into the gapless CFL 
phaseJ23|- If, instead, we choose a larger Aq, as in 
Ref. I3, we find that our results become markedly 
more A-dependent, which is a good diagnostic for 
model-dependence. 

4. We work to leading nontrivial order in Ai, A2, A3, 
/ie, /i3 and fig. This should be a good approxima- 
tion, as all these quantities are small compared to 
/i. 

5. We neglect the anti-particles. This simplifies the 
numerics by discarding physically unimportant de- 
grees of freedom, but one must be cautious with 
this truncation. It introduces cutoff-dependent 
terms in our free energy, including some that de- 



pend on the chemical potential and therefore in- 
troduce cutoff-dependence in the corresponding 
charges. For our purposes this is not important, 
firstly because we always present free energy dif- 
ferences relative to neutral unpaired quark mat- 
ter, and secondly because we only care about elec- 
tric and color charges that have zero trace over 
all fermion species, and for these the cutoff depen- 
dence cancels out. However, a non-traceless charge 
like baryon number would have an incorrect cutoff- 
dependent value when calculated in this approxi- 
mation. 
6. We ignore meson condensation in both the CFL 
and gCFL phases. 



We expect that these approximations have quantitative 
effects, but none preclude a qualitative understanding of 
the new phase we shall describe. 

We now give the explicit expression for detS*"^, after 
having implemented the approximations above. As de- 
scribed in section rTBl we use a color- flavor basis in which 
the gap matrix |(2Jl is conveniently block-diagonal. Since 
the chemical potential and mass are diagonal in color 
and flavor, the full inverse propagator (|16|l is then also 
block-diagonal in color-flavor space. This means we can 
break the determinant in equation H21|) into four more 
manageable pieces: 



detS~^{po,p) = 



-"-^ru ^ rd ^rs gs 

gd gu bu bd 
bs 



(23) 



We find that the 2x2 determinants are 



Drd = IQfJ.'^iiprd - P - iP()){f^gu - P + iPo) + ^l) 
{(P'rd-P + iPo){P'gu-p-m)+ Aj), 



D^ 



bu 



D„ 



bd 



16/i**((/irs - P - ipo){lJ-bu ~P + ipo) + A^) 
((/i^s -P + ipo)(/ifcn -p- ipo) + A|), 

16/i**((/i6d - p - «Po)(/*gs -P + ipo) + A^) 
{{P'bd - P + ipQ){Pgs -p~ipo) + Af), 



and the 3x3 determinant is 



(24) 



I?™ = a(2/.)6( b{de-Al)icf-Al)-cdfAl 



gd 
b 



-dAfAl - defAj + fAlAl 



-(2/.)6(^ AlideAl~4AlAl + efAl) 
+c{dAi + fAlAl) 
+b{eAlAl + c{AlAl - deAj - e/Ai)) 



(25) 



(26) 



where for compactness we assign 

a = ^J-ru-p + ipo 

b = - i-i.ru +p + ipa 

c = iJ-gd-p + ipo 

d = -/igd +P + iPo 

e = iibs-p + ipo 

f = -j-ihs +p + ipo 

and where we have dropped the superscript on the "ef- 
fcctive quark chemical potentials" , given by 

Uru = M ^ 3^6 + 2^3 + 3^8: 
l^ld = M + 3^6 - 2^3 + 3^8, 
/^bf - f^+^t^e- |M8 - M2/(2A^), 

Mrd = Ai + iMe + ^M3 + 5^8, 

Mg« = M - 3^6 - 2/^3 + 3/^8, ,2Y) 

iif, = /i + i^e + 5^3 + 3M8 - Af.V(2A'), 



cff 

^^bu 


= ^- l^e- |m8, 






= /i + i/Xe - i/X3 + i/Z8 - 


- M^/{2fi) 


,,eff 


= /i + i^e - 1^8 • 





These expressions explicitly show how we treat the 
strange quark mass as a shift in the chemical potential 
of the strange quarks. In evaluating these determinants, 
we have extensively used the identity 



det 



A B 
C D 



det{A)det{D~CA-^B) 



for the determinant of a block matrix. 

The numerical task is now explicit. We find the quasi- 
particle dispersion relations e{p) by finding the zeros of 
(|23|l . viewed as a polynomial in pq. We then perform the 
integral in (|20|l numerically, and obtain fl. We evaluate 
the partial derivatives of the free energy required in the 
neutrality conditions ||SJ) and the gap equations H22|) nu- 
merically as finite differences, with differences 0.1 MeV 
in the relevant chemical potential or gap parameter. 

As a check, we have also done the calculation of fi and 
its partial derivatives by evaluating both the p and po 
integrals numerically, never writing the latter as a Mat- 
subara sum. In this alternative calculation, we were able 
to evaluate the partial derivatives in and 122|) analyt- 
ically. 

The n that we obtain is cutoff dependent, but its par- 
tial derivatives ((SJ and (|22|l are insensitive to variations 
in A in the sense described above, namely as long as the 
coupling is changed to keep Aq fixed upon variation in 
A, and as long as A is kept well above /x. Furthermore, 
we are only ever interested in free energy differences be- 
tween phases. When we evaluate the differences between 



the fi for unpaired, (g)CFL, and (g)2SC quark matter, 
we find that all such free energy differences are insensitive 
to the cutoff, as they should be since these differences all 
reflect physics near the Fermi surfaces. Because we are 
only interested in free energy differences, in evaluating fi 
we make the numerical integral better behaved by sub- 
tracting the appropriate expression for neutral unpaired 
quark matter within the integrand. 

The solutions of the system of gap and neutrality equa- 
tions depend on three parameters: fi, Ms and Aq. Our 
purpose is to understand the effect of Mg on CFL pair- 
ing, and these effects are controlled by the relative size 
of M"^ / PL and the gap parameters A^, whose overall mag- 
nitude is set by Aq. It is therefore better to think of the 
three parameters in the problem as /i, M^//i and Aq. In 
compact stars, ^ increases and Ms presumably decreases, 
meaning that M^//j, decreases as one approaches the cen- 
ter of the star. For simplicity, we set the overall energy 
scale in our calculation by fixing /i = 500 MeV, which is 
reasonable for the center of a neutron star, and vary Mg 
in order to vary Mg/fi. We have confirmed that as long as 
we choose a Aq that is small enough that the transition 
© occurs where M"^/ ^j? corrections are under control, 
this transition occurs very close to M^//j, « 2A, where 
A is gap parameter on the CFL side of the transition. 
The authors of Ref. [3 have confirmed that this result 
continues to be valid even for Aq as large as 100 MeV, 
where the approximations are not as well controlled. We 
quote results only for Aq — 25 MeV, which is within the 
plausible range of values that Aq may take in nature [l| 
and for which our calculation is clearly under control. 
Although we have obtained our results by varying Mg at 
fixed /i, we typically quote results in terms of the impor- 
tant combination M'^/ii. 



III. RESULTS 
A. Domain where gCFL is favored 

In Figs, n and |21 we show the gap parameters and 
chemical potentials as a function of Af^//i,, for Ag = 
25 MeV. Fig. O shows the free energy. We see a con- 
tinuous phase transition occurring at a critical M^ that, 
in our model calculation with /i = 500 MeV, lies be- 
tween Ms = 153 MeV and Ms = 154 MeV, i.e. at 
{M"^ / li)c ~ 47.1 MeV. This agrees exceedingly well with 
the expected value 2 A from Eq. lO, since on the CFL 
side of the transition Ai = A2 = A3 = 23.5 MeV. For 
Mg/fi < {M^/fj,)c, the CFL phase is favored, with all 
three gaps equal to each other within our approxima- 
tions. If we improve upon our approximate treatment 
of Ms, we expect Ai = A2 with these gap parameters 
slightly smaller than A3, because Ai and A2 describe 
pairing between quarks with differing Fermi velocities, an 
effect of Ms that we are neglecting because it is known 
to be small in other contexts |2g. (Indeed, it proves to 
be a few percent effect also in the present context |23i-) 
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FIG. 1: Gap parameters A3, A2, and Ai as a function of 
Ms/fi for fj. = 500 MeV, in a model where Ao = 25 MeV (see 
text). At M^/fi « 47.1 MeV (vertical dotted line) there is 
a continuous phase transition between the CFL phase and a 
phase that we shall identify below as the gapless CFL phase. 
We find gapless CFL phase solutions up to Ms/fi « 144 MeV. 
But, we shall see in Fig. ^ that above M^/fi « 130 MeV 
(which we denote here with a vertical dash-dotted line) un- 
paired quark matter has a lower free energy than the gapless 
CFL phase. 





10 


r—i 




> 














^H 


■10 


01 




U 




a 


■20 


U 




01 











■30 


>-, 




sr 






■40 


a 




M 





-50. 




50 75 

MJ[i [MeV] 



150 



FIG. 3: Free energy of the CFL/gCFL phase, relative 
to that of neutral noninteracting quark matter and that of 
the 2SC/g2SC and 2SCus phases, discussed in Section [ill El 
There is a CFL^gCFL transition at {M'^/fj.)c ^ 47.1 MeV, 
(vertical dotted line), at which the free energy and its slope 
are continuous, indicating that the transition is not first or- 
der. If we neglect the possibility of other phases (for example 
g2SC) we would conclude from this figure that there is a first 
order transition gCFL^unpaired at Ms /fj. ~ 130 MeV (ver- 
tical dash-dotted line). 



-M^/(2/i) and the free energy is 
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FIG. 2: Chemical potentials ^e, fJ-3 and ^s as a function of 
Mg /jj. in the CFL/gCFL phase for the same parameters as in 
Fig. The effects of electrons on the free energy have been 
included in the calculation, as will be discussed in more detail 
below. We see that the gapless CFL phase has jj,e > 0, mean- 
ing that it has a nonzero density of electrons. Perhaps the 
most physically relevant order parameter for the CFL/gCFL 



phase transition is the electron number density n^ 



'Me 



For M^/^ < {Mg/fi)j~our results agree with the small- 
Ms expansion of Ref. [S|, where /^a = /ie = and /xg = 



^ ^ neutral 
CFL 



3// 



47r2 

3A2/^2 



l-121og(Af,/2M) 4 
327r2 



^ ^ neutral 

unpaired 



3Mf - 48A2^2 
167r2 



(28) 



As the density decreases (i.e. as M'^/n increases) 
through the CFL^gCFL transition, the gap parameters 
split apart, with A3 increasing slightly and A2 and Ai 
dropping significantly, with Ai dropping faster than A2. 

We have verified that M'^/fj.A is the relevant di- 
mensionless quantity by changing the coupling strength, 
i.e. picking a different Aq (gap at Mg = 0). The critical 
point {M^/fi)c changes as predicted by ®. Furthermore 
we checked the robustness of our results upon variation 
of the cutoff A, observing changes of only a few percent 
in the value of M'^ / fiA at the transition upon changing 
A by up to 20% while keeping Aq fixed. 

Fig. O confirms that the slope of the free energy is 
continuous at the CFL/gCFL transition, indicating that 
it is not first order. We have not determined the or- 
der of the transition, because evaluating higher deriva- 
tives of the free energy with respect to M"^ / fi is not nu- 
merically feasible. The most physically relevant order 
parameter is the electron density rif. ~ /Xg, which is of 
course equal in magnitude to the electric charge density 
of the quarks. This increases above the transition like 
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FIG. 4: Dispersion relations at Mf//i = 80 MeV (with /i 
and Ao as in previous figures) for the quasiparticles that are 
hnear combinations of gs and bd quarks (dashed hnes) and 
for the quasiparticles that are linear combinations of bu and 
rs quarks (solid lines). There are gapless gs-bd modes at 
pY = 469.8 MeV and p^ 509.5 MeV, which are the bound- 
aries of the "blocking" |29l 13(1 or "breached pairing" |31ll re- 
gion wherein there are unpaired bd quarks and no gs quarks. 
One bu-rs mode is gapless at p = 475.6 MeV with an almost 
exactly quadratic dispersion relation that we shall discuss be- 
low. 



Ue ~ [(Af^//i) — [M"^ / fi)^ , suggesting a fourth order 
phase transition! This argument neglects the small elec- 
tron mass and, furthermore, it neglects the fact that, as 
we shall see in Eq. (|33|) . there is also a nonzero num- 
ber density of neutral unpaired quark quasiparticles that 

1 /9 

grows like [{M"^ / ^) - (M^^//i)c] . Although because 
these unpaired quasiparticles are neutral they are less 
important phenomenologically, this does suggest that the 
transition is second order, as in the analysis of Ref. ^28] . 

If we had used a simpler ansatz in which the gap pa- 
rameters were constrained to one value A = Ai = A2 = 
A3, then the CFL phase would have remained artificially 
stable above the critical value of M^/fi. From Eq. (|28|l . 
its free energy would be higher than that of gCFL, rising 
to equality with that of unpaired quark matter at a value 
of M^/fi around 90 MeV. (The precise value depends on 
how the common A changes with Ms.) 

Of course, we actually used the more general ansatz 
© that allows the A's to differ. We found that the CFL 
phase becomes unstable and is replaced by the gCFL 
phase, in which the gaps have very different values, so 
the simplified analysis of Ref. [g does not apply. The 
free energy of the gCFL phase crosses that of unpaired 
quark matter at M^/n « 130 MeV. This phase transition 
is first order, and we are able to follow the metastable 
gCFL phase up to M^ //x = 144 MeV where, as we shall 
explain below, it ceases to be a solution. 
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FIG. 5: Dispersion relations at M^/^i = 80 MeV for the 
two quasiparticles that are linear combinations of rd and gu 
quarks (dashed and dotted), and for the three quasiparticles 
that are linear combinations of ru, gd and bs quarks (solid). 
These five quark quasiparticles all have gaps throughout the 
CFL and gCFL phases. 



B. The nature of the gCFL phase 

Up to this point we have not justified our use of the 
name "gapless CFL" for the new phase that replaces 
the CFL phase at M^/fJ, > 2A. We have given model- 
independent arguments to expect that it will contain un- 
paired bd quarks, but now we describe its properties in 
more detail. In calculating the free energy (|20|l of the 
Cooper-paired quark matter we automatically obtain the 
quasiquark dispersion relations (|21|l . so we can see what 
gapless modes exist. These modes are important because, 
at the temperatures T < keV characteristic of neutron 
stars, only the lightest modes will contribute to trans- 
port properties. 

In Fig.^we show the dispersion relations for the rs-bu 
and gs-bd 2x2 blocks in the quasiquark propagator, at 
Mg /fJ, = 80 MeV. We see immediately that there are 
gapless modes in both blocks, justifying our name for 
this phase. Before moving on to a detailed discussion 
of the physical properties of the gCFL phase, we should 
note that the phenomenon of gapless superconductivity 
is well known, at least theoretically. It was first sug- 
gested by Sarma |32l | who worked in a context much like 
our gs-bd block in isolation, and found that the gapless 
superconducting phase is never stable. Alford, Berges 
and Rajagopal found a metastable gapless color super- 
conducting phase in Ref. |33j| , but this phase was neither 
electrically nor color neutral. The key observation was 
made by Shovkovy and Huang ^^, who discovered that 
when the constraints of electric and color neutrality are 
imposed on the 2SC phase in two-flavor QCD, there are 
regions of parameter space where a gapless color super- 
conducting phase is stable. Following their nomenclature 
(they described a "gapless 2SC phase") we refer to the 
phase that we find above M^ as the "gapless color-flavor 
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locked phase" [3. 

Gapless two-flavor color superconductivity was also 
studied in Ref. [311, building upon prior work done in 
a cold atom context [SJI- These authors analyzed pair- 
ing between a heavy and a light quark, akin to gs and 
bd, in the case in which the gs quarks are nonrelativistic. 
They find that a gapless phase (they describe the block- 
ing region as a region in which pairing is "breached" ) is 
stable if the relative density of the two species is held 
fixed. 

Note that in our three-flavor calculation, both the gap 
equations (|22|) and the neutrality conditions ||HJ) couple 
all nine quarks. Although the single particle dispersion 
relations can be analyzed for the gs and bd quarks in 
isolation, and are qualitatively similar to those obtained 
in Refs. |2lL l31| in a two-flavor setting, the implications 
of neutrality are more subtle in our three-flavor context 
as we shall explain below. 

Each of the dispersion relations in Figs. 0] and [S] de- 
scribes an excitation with well-deflned Q, although the 
sign of Q changes at momenta where the dispersion re- 
lation is gapless. Beginning with an example with no 
gap, the upper solid curve in Fig. ^describes excitations 
that are linear combinations of rs particles and bu holes, 
both with Q — —1. The lower dashed curve in Fig.Qlhas 
clearly visible momenta pj'' and P2'' where it is gapless, so 
we use this as an example of "sign change" even though 
it describes Q = quasiparticles: to the left of pY, it 
describes (7s-holes with a very small admixture of bd par- 
ticles; to the right of P2'', it describes bd particles with a 
very small admixture of gs holes; but, between p^"^ and 
^2'^ it describes excitations that are superpositions of bd 
holes and gs particles. 

In the CFL phase, once we take into account the 
explicit symmetry breaking introduced by the strange 
quark mass and elcctromagnctism, the unbroken sym- 
metry is reduced from the diagonal SU(3)l+r+c to 
U{1)q X C/(l) _12J. The last C/(l) corresponds to "color 
-|- flavor hypercharge" and may be spontaneously broken 
by meson condensation _13j . The gapless CFL phase has 
the same symmetry as the CFL phase, and it will there- 
fore be interesting to investigate the possibility of meson 
condensation in the gCFL phase. The effective theory for 
the Goldstone bosons alone will have the same form as 
in the CFL phase, albeit with new contributions to their 
masses coming from the differences between the values of 
the three A|. And, furthermore, the gapless quasiparti- 
cles must be included in the low energy effective theory. 



1. Dispersion relations, gapless modes, and neutrality 

As will soon become clear, the 3x3 block in the pairing 
pattern Q plays a minor role: its quasiparticles are al- 
ways gapped, so we mainly discuss the three 2x2 blocks. 
In general, when two species of massless quarks undergo 
s-wave pairing with gap parameter A, the dispersion re- 



quark pair 



SUcH 



SfJ-clt in 
electronless CFL 



rd-gu \{^ie+^z) f^e 

rs-bu i(^e + i^3 -f /i8 - ^Ms/m) Me - |Afs /m 



gs-bd 



;(iM3-M8 + |M,V/i) 



Mf/2/. 



TABLE I: Chemical potential splittings for the 2x2 pairing 
blocks. (Sficft and /i, which is not tabulated, are defined in 
each row such that the efi^ective chemical potentials of the two 
quarks that pair are jl±SfieS-) The middle column gives (5/ioff 
for general values of the chemical potentials fj.e, ^3 and fj,s. 
In the last column, it is understood that as fi^ is varied, ^3 
and /is "follow it" in such a way that varying fie corresponds 
to varying /Xq, tracking degenerate Q-neutral solutions for 
electron- less CFL quark matter. 



lations of the two resulting quasiparticles are 



E{p) 



Sfi±^ip~fl)^ + A^ 



(29) 



where the individual chemical potentials of the quarks 
are p, ± Sfj,. As long as the chemical potentials pulling 
the two species apart are not too strong. Cooper pairing 
occurs at all momenta: 



pairing criterion: \Sfi\ < A 



(30) 



However when this condition is violated there are gapless 
{E = 0) modes at momenta 



Pgaplcss 



= /2 ± V V - A2 



(31) 



and there is no pairing in the "blocking" or "breached 
pairing" region between these momenta |2j |23, yfl IsJ 
I2JI . (The identification of the boundaries of a blocking 
region with locations in momentum space where a disper- 
sion relation is gapless is discussed with considerable care 
in Ref. [30j, which considers a more complicated setting 
in which rotational symmetry is spontaneously broken 
and the blocking regions are not spherically symmetric. 
Such blocking regions were analyzed previously in Refs. 
^9"!.) The pairing criterion H30() can be interpreted as 
saying that the free energy cost 2A of breaking a Cooper 
pair of two quarks a and b is greater than the free energy 
25/i gained by emptying the a state and filling the b state 
(assuming that 6fi pushes the energy of the a quark up 
and the b quark down) '25'. In the blocking region, we 
find unpaired 6 quarks and no a quarks. 

We wish now to apply these ideas to the 2x2 pair- 
ing blocks in three-fiavor quark matter, first in the CFL 
phase. As described above, neutrality is imposed via 
chemical potentials fie, ^J■3, ^J■SJ ^^nd in the CFL phase 
the leading effect of the strange quark mass is an ad- 
ditional effective chemical potential — Af^/2/i for the 
strange quarks. The splittings of the various pairs are 
then as given in the middle column of table HJ 

Electrons will play a crucial role in understanding the 
gCFL phase, but it is fruitful initially to consider matter 



11 



40 



a> 



20 



^ 



-20 



I ' ' ' ' I 



rd-gu unpairing line 




gs-bd unpairing line 



25 



50 



75 
M:/^ [MeV] 



100 



125 



150 



FIG. 6: Unpairing lines for the same parameters as used 
in Fig. If electrons are neglected, then the upper and 
lower curves bound the region of fie where neutral solutions 
to the gap equations are found. These solutions are all Q- 
insulators. Taking electrons into account, the correct solution 
is the dashed line: in the CFL phase jj,e — 0, while the gCFL 
phase corresponds to values of fie below but very close to the 
rs-bu unpairing line. gCFL is a Q-conductor both because 
of the nonzero electron density and because of the ungapped 
Q-charged rs-bu quasiparticles. 



consisting only of quarks, which wc can do by sending 
the electron mass to infinity. In the absence of electrons, 
at each M'^/fi, there is a plateau in the free energy of 
neutral CFL (or gCFL) solutions: if we vary the chemical 
potential that couples to Q charge, 



^^Q 



I (Me 



M3 



^Ms) 



(32) 



while keeping constant the gap parameters A^ and the 
two orthogonal combinations of chemical potentials, then 
over a range of fiQ the free energy does not change and 
we have a family of neutral stable solutions to the gap 
equations. This indicates that, in the absence of elec- 
trons, both the CFL and gCFL phases are Q-insulators. 
On this plateau, all Q-charged quasiparticles remain 
gapped: (|Sn|l is obeyed for the {rd,gu) and {rs, bu) 2x2 
quark pairing blocks. At the edges of the plateau, some 
Q-charged quasiparticles become gapless, the material 
ceases to be Q-neutral, dQ./dfiQ ^ 0, and the free energy 
is no longer independent of changes in /iA . The range of 
/ig that defines the plateau is therefore the band gap for 
the CFL/gCFL insulator. In Fig. we show the unpair- 
ing lines for each 2x2 quark pairing block. The rd-gu line 
and the rs-bu line bound the plateau region. Although 
the vertical axis is labelled "fie" , it actually corresponds 
to variation in /j,a, since we varied {fie, fis, fis) by a mul- 
tiple of (1, 1, i). In the CFL phase, this corresponds to 
keeping fj,^ = fie and fig = \{^e — M^/fi) while varying 

^J■e■ 

We see that (g)CFL matter exists in a wedge, between 
the rd-gu unpairing line and the rs-bu unpairing line. 
From table ^ we can see that the bd- gs-nnpairing line 



is vertical because the bd and gs quasiparticles are Q- 
neutral, so their splitting depends only on M'^/fi and 
not on fie- This unpairing line has a different character 
than the other two. Rather than bounding the band-gap 
within which solutions are found, it separates the CFL 
and gCFL phases. CFL is stable only up to a critical 
value of M'^/fi, where the the gs-bd pairs break. 

At the lower (rs-6u-unpairing) line, fi^ is large enough 

that the bu and rs quarks, which have Q — -\-l and 
Q = — 1 respectively and which pair with gap param- 
eter A2, no longer pair completely: it is energetically 
favorable to create a new blocking region of unpaired bu 
quarks. At this Q-electrostatic potential, the CFL Q- 
insulator breaks down, unpaired bu quarks with Q — -f 1 
are created, the free energy is no longer /ig-independent, 
and in fact the neutrality conditions and gap equations 
are no longer satisfied. 

At the upper (rd-(7u-unpairing) line, fip is so low that 

the rd and gu quarks, which have Q = —1 and Q ~ +1 
respectively and which pair with gap parameter A3, no 
longer pair completely, and it is energetically favorable 
to create a new blocking region of unpaired rd quarks, 
and once again no solution is found. 

At M^/fi = 143 MeV, which is so large that the gapless 
CFL phase is anyway already metastable with respect to 
unpaired quark matter, the two boundaries cross, mean- 
ing that no gapless CFL solution can be found. 

So, in the absence of electrons, we can find stable so- 
lutions of the gap and neutrality equations everywhere 
between the rs-bu and rd-gu curves in Fig. To the 
left of the gs-bd unpairing line this is the CFL phase, a 
Q-insulator with no gapless quasiquark modes. To the 
right of that line we have the gCFL phase, again a Q- 
insulator, in which all Q-charged modes are gapped, but 
there are Q = gapless quasiparticles. 

We now restore the electrons, setting their mass to 
zero. In the CFL region, the system is forced to fie = 
(dashed line in Fig. EJ I23. However, at the transition 
point to gCFL, where the gs-bd pairs break, we find that 
the neutrality requirement forces us over the line where 
rs-bu pairs also begin to break. The result is that as 
M^ / fi increases further, the system maintains neutrality 
by staying close to the rs-&M-unpairing line, where there 
is a narrow blocking region in which there are unpaired 
bu quarks. Their charge is cancelled by a small density 
of electrons. We analyze this quantitatively below. 

We see that real-world gCFL quark matter is a con- 
ductor of Q charge, since it has gapless Q-charged quark 
modes, as well as electrons. The rd and gu quarks, which 
are insensitive to the strange quark mass, remain robustly 
paired and the Q-neutral bd and gs quarks develop a large 
blocking region as the system moves far beyond their un- 
pairing line. The neutrality requirement naturally keeps 
the system close to the rs-6u-unpairing line, following 
the dashed line in Fig. [^ so these quarks have a very 
narrow blocking region and an almost quadratic disper- 
sion relation (see below). Although U{1)q is unbroken in 
the gapless CFL phase, the presence of electrons and un- 
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paired bu quarks makes this phase a Q-conductor. This 
is in contrast to the CFL phase, which is a Q-insulator 
with no gapless quasiquarks and no electrons. 

The gapless quark quasiparticles occur in the gs-bd and 
rs-bu sectors. Since these will have a dramatic effect 
on transport properties, we now discuss them in greater 
depth. 



phase transition, the long wavelength physics at the crit- 
ical point is 4-dimensional rather than 3-dimensional as 
at a finite temperature transition. 



The rs-bu sector 



2. The gs-bd sector 

In a typical part of the gCFL phase space, the Q- 
neutral gs-bd sector is well past its unpairing line, and 
there is a large blocking region between momenta pj'^ 
and pj'^ at which there are gapless excitations, as shown 
in Fig. 21 In the blocking region p^"^ < p < p^2 there are 
bd quarks but no gs quarks, and thus no pairing. We 
have confirmed this by direct evaluation of the difference 
between the number density of bd and gs quarks, show- 
ing this to be equal to the volume of the blocking region 
in momentum space. 

Note that even though there is no pairing in the ground 
state in the blocking region, the dispersion relations are 
not trivial. Because the states obtained via the two dif- 
ferent single particle excitations that are possible (adding 
a gs quark or removing a bd quark) mix via the Ai con- 
densate, the two dispersion relations exhibit an "avoided 
crossing" between p^-^ and p^-^ . If we neglect the mixing 
among the excitations introduced by Ai, the gapless ex- 
citations just above (below) pl^^ ^-re bd quarks (holes) and 
those just above (below) pj'' are gs quarks (holes). 

It may seem coincidental that the value of M"^ / jjl at 
which the CFL phase becomes gapless is the same as the 
value at which Ai and A2 separate in Fig. ^ Although 
we do not see a profound reason for this, it is certainly 
not a coincidence. The CFL^gCFL transition is trig- 
gered by the instability of the CFL phase that occurs 
when a gs-bd quasiparticle dispersion relation goes gap- 
less, indicating the instability towards gs-bd unpairing 
and the opening up of a blocking region in momentum 
space, filled with unpaired bd quarks and with no gs- 
bd pairing. Consequently, one of the terms in the Ai 
gap equation — that corresponding to the gs-bd block 
— is reduced in magnitude because its integrand van- 
ishes within the blocking region. This reduction in the 
support of the Ai gap equation integrand causes Ai to 
drop. 

The "thickness" of the bd blocking region can be con- 
sidered an order parameter for gCFL: for M|//i below 
the critical value there is no blocking region. Just above 
the critical value we can use the results of table Q] and 
(|31|l to show that 



J^d _ bd 
P2 Pi 






typical behavior for a second order phase transition. 
Because we are analyzing a zero temperature quantum 



As discussed above, the gCFL phase remains neu- 
tral by crossing the rs-bu unpairing line, and developing 
enough unpaired bu quarks to cancel the Q charge of the 
electrons. The electrons contribute (— /if/127r^) to the 
free energy, so Q-neutrality can be maintained as long as 
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{pTf ~ JPl) 
37r2 



(34) 



where p\^ and pj" bound the blocking region of unpaired 
bu quarks. The condition (|34|l implies that 



^2- 



pT) 



3p2 



(35) 



where p is the average of pj" and ^2"- ^t M^/^ = 
80 MeV, where fie — 14.6 MeV at the lower curve in 
Fig. El this imphes (p^" -pj") = 0.0046 MeV! Indeed, in 
Fig-Elthe separation between p\^ and P2" is invisible, and 
the dispersion relation appears to be quadratic about a 
single gapless point. To resolve the separation between 
p\" and P2"; '^^ did calculations assuming 200 and 500 
"flavors" of massless electrons. In these cases, (p2"~Pi") 
~ 1 MeV and ~ 3 MeV, in very good agreement with the 
above argument. Returning to our world with its single 
electron species, because (pa" ~Pi") is so small, the value 
of /ig at the true Q-neutral solution is very close to that 
given by the lower curve in Fig. And, the gaps are very 
close to those found in a calculation done in the absence 
of electrons. 

From Eq. (|29|l . the maximum in the quasiparticle en- 
ergy between the two gapless momenta is i^max — l^fJ- ~ 
A I , so from (|31|l we can express this in terms of the width 
of the blocking region: 4(|(5^|-|- A)i?„iax = (P2— Pi)^- For 
the rs-bu quarks, the blocking region is always very nar- 
row, so £;max ~ (p'r - Pi")V(8A2) which from (jSHJ is 
a small fraction of an electron volt at M'^ /fj, = 80 MeV. 
Thus, at any astrophysically relevant temperature, the 
rs-bu dispersion relation can be treated as quadratic 
about a single momentum at which it is gapless. In- 
deed, even at Af^ //i = 130 MeV where the gapless CFL 
phase ceases to be the ground state, /ig = 40.3 MeV, 
(P2" — p5") ^ 0.1 MeV and the peak of the dispersion 
relation between pj" and pj" is at about 50 eV. The re- 
quirement of Q-neutrality naturally forces this dispersion 
relation to be very close to quadratic, without requiring 
fine tuning to a critical point. 
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4- The rd-gu sector and the 3x3 ru-gd-bs block 

In Fig. [SI we show the dispersion relations for the quasi- 
particles in the rd-gu sector. One of these becomes gap- 
less at the upper boundary of the wedge in Fig.EJ but we 
have seen that in the presence of electrons, the neutral 
gapless CFL solution is never near this upper boundary. 
Therefore, these dispersion relations are always gapped, 
as in the figure. In Fig. [S] we also show the dispersion re- 
lations for the three quasiparticles from the 3x3 block. 
These quasiparticles carry zero Q charge and they always 
have nonzero gap. Their smallest gap becomes very small 
near the rightmost tip of the gCFL wedge region in Fig.EJ 
but is always greater than 1 MeV in the region in which 
gCFL is favored. 



C. The gCFL EVee Energy Function 

In the previous subsection, wc have used the dispersion 
relations to delineate the unpairing lines which bound 
the ranges in ^q where, in the absence of electrons, Q- 
insulator solutions are to be found and which separate 
the CFL and gCFL phases. Here, we sketch the behavior 
of the free energy fl in the vicinity of solutions to the 
gap and neutrality equations, and see how this behavior 
changes at the unpairing lines. 

In Fig. [3 we study the free energy in the vicinity of 
a gapless CFL solution not far above the CFL— >gCFL 
transition. We have neglected electrons in making this 
plot; the change from including them would be invisi- 
ble on the scale of the plot. We plot the free energy 
upon variation of gap parameters while keeping /x's fixed 
(dashed curves in Fig. EI), and we also plot the "neutral 
free energy" (solid curves) obtained by varying gap pa- 
rameters about the solution while solving the neutrality 
conditions anew for each value of the gap parameters. 
We see that the solution is a minimum of the neutral 
free energy, confirming that we have succeeded in finding 
a stable neutral solution. However, the solution is not at 
a minimum of the free energy upon variation of the gap 
parameters while keeping /x's fixed. 

We see in the top panel of Fig. [7| that the solution is 
found at a local maximum of the dashed curve describ- 
ing variation of Ai at fixed /x's. Shovkovy and Huang 
described similar behavior in the gapless 2SC case in 
Refs. |2^[23, and suggested that this is a characteris- 
tic of gapless superconductivity. We find this not to be 
the case: deep in the gCFL phase, at M^/ fi = 80 MeV 
rather than the M^/n = 51.2 MeV of Fig. [3 wc find that 
the gCFL solution is a local minimum of both the dashed 
and the solid curves in the analogue of the top panel of 
Fig. 13 That is, we find an onset of the behavior seen 
in the top panel of Fig. [3 as we cross the CFL— >gCFL 
transition, namely the gs-bd unpairing line: as a gs-bd 
blocking region begins to open up, the solution goes from 
being a local minimum of the dashed curve to being a lo- 
cal maximum. However, we find that the dashed curve 
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FIG. 7; These figures show the free energy Q in the vicinity 
of the gapless CFL solution for M^/fj, = 51.2 MeV. In each 
panel, the dashed curve is obtained by varying one of the 
gap parameters (Ai in the top panel, A2 in the middle; A3 
in the bottom) while keeping the other two gap parameters 
and the chemical potentials /^e, /is and fig fixed. The free 
energies are measured relative to that of the solution. The 
solid curve in each panel depicts the "neutral free energy", 
obtained by varying one gap parameter, keeping the other 
gap parameters fixed, and solving the neutrality conditions 
anew for each point on the solid curve. 



does not persist in this shape as the gs-bd blocking region 
expands. The onset of gaplessness is characterized by a 
dashed curve as in the top panel of Fig. [3 but gaplessness 
itself need not be. 

In the middle panel of Fig. [3 we find that the solution 
is at a point of inflection with respect to variation of 
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FIG. 8: Same as the middle panel of Fig. |7| except that /i^ 
has been increased by 2 MeV while changing /13 and /ig so 
as to make this a shift in jIq. This means that, neglecting 
electrons, we are now exploring the change of the free energy 
and the neutral free energy upon variation of A2 about a 
gCFL solution that is in the interior of the wedge in Fig |S| 
rather than at its bottom boundary. 



A2 at fixed ^'s. We find that the gCFL solutions at all 
values of Mlj \i above [M'^j ij)c are at points of inflection 
of this sort. This arises because a gCFL solution is forced 
by the neutrality constraint to be very close to the bu- 
rs unpairing line. In Fig. [S] we replot the middle panel 
of Fig. [7| after increasing /Ze by 2 MeV while varying /is 
and /xg so that only /iA changes. This means that we 
have taken a 2 MeV step upwards in Fig. |H| away from 
the bu-rs unpairing line. And, we see that the solution 
is now a minimum with respect to variation of A2 at 
fixed /x's. The point of inflection has resolved itself into 
a minimum and a maximum, with the solution at the 
minimum. Thus, the point of inflection in the dashed 
curve does indeed occur at the bu-rs unpairing line. 

Note that at M"^ / ji = 80 MeV, once we have taken an 
upward step away from the bu-rs unpairing line in Fig. (BJ 
obtaining the analogue of Fig. |H1 the gapless CFL phase 
solution is now a local minimum of both the dashed and 
solid curves for variation in the Ai, A2 and A3 directions. 

If we take a step in the "wrong direction" in Fig. El 
downwards from the bu-rs unpairing line, the point of 
inflection in the middle panel of Fig. |7| vanishes and the 
dashed curve becomes monotonically increasing, indicat- 
ing that there is no solution to the A2 gap equation to 
be found at these values of the /i's. In the presence of 
electrons, the neutrality conditions are satisfied just be- 
low the bu-rs unpairing line in Fig. |^ and the depen- 
dence of the free energy on A2 is slightly modified so 
that the point of inflection in the middle panel of Fig. [7| 
occurs where the neutrality conditions are satisfied. We 
have confirmed this in calculations done with 200 and 500 
species of electrons; with a single species as in the real 
world, the changes in Fig. |7|are invisible on the scales of 
the plot. 

Finally, with respect to variation of A3 , the solution is 
a local minimum of the dashed curve in the lower panel 
of Fig. \7\ However, we have verified that if we move 




FIG. 9: Schematic illustration of conditions for the occur- 
rence of mixed phases. Free energy Q for two phases A and B 
is shown as a function of some chemical potential /ii . Charge 
Qi — —dQ/dfii is given by the slope. Squares mark the neu- 
tral points. Panel (a): at the neutral point for each phase, the 
other phase has lower free energy, so there is a point (black 
dot) where the two phases can coexist with the same pressure 
and opposite charge, with lower free energy than either neu- 
tral phase. Depending on Coulomb and surface energy costs, 
a mixed phase may exist there. Panel (b): phase B has higher 
free energy than phase A at the point where A is neutral. At 
no point do the two phases coexist with opposite charge, so 
no mixed phase is possible. 



sufficiently upwards in Fig. as to run into the rd-gu 
unpairing line, then the dashed curve in the lower panel 
exhibits a point of infiection (while that in the middle 
panel has a robust minimum.) 



D. Mixed Phase Alternatives 

Up to this point, the phases we have discussed have 
been locally neutral with respect to all gauge charges. 
However, it is well known that neutrality can also be 
achieved in an averaged sense, by charge separation into a 
mixture of two oppositely charged phases. This is shown 
schematically in Fig. |^ which shows generic free energy 
curves ft{iJ,i) for two phases A and B. The free energy 
must be a concave function of the chemical potential, 
since increasing /Zi increases the charge Qi = —dVl/d^i. 
There are then two possible situations. In one (Fig.jSJD) 
there is no coexistence point and hence no mixed phase 
is possible. In the other (Fig. |5Ji) there is a coexistence 
point of oppositely-charged phases, and its free energy is 
lower than that of either neutral phase, so if Coulomb 
and surface energy costs are low enough then a neutral 
mixed phase will be free-energetically preferred over ei- 
ther homogeneous neutral phase. 

We now consider possible gCFL-|-unpaired mixed 
phases. (Note that for (M^/fi) > {M^/^i)c a 
CFL-|-unpaired mixture is not possible because there is 
no CFL solution, charged or neutral.) For the unpaired 
and gCFL phases, the free energies are of the form shown 
in Fig. |5^. At the values of {nei ^^^,|^s) that make one 
phase neutral, the other phase has lower free energy. 
Thus there is a value of {ii^, fJ-s, fJ-s) "between" that for 
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neutral gCFL and that for neutral unpaired quark mat- 
ter, where the two can coexist with opposite color and 
electric charge density. However, a mixed phase is not 
favored in this case because each component would have 
net color charge, and color is a gauge symmetry with a 
strong coupling constant, so this mixed phase would pay 
a huge price in color-Coulomb energy. (For similar ar- 
guments applied to systems with no gauge symmetries, 
where the initial conclusion that a mixed phase is favored 
is the correct one, see Refs. [35|,|36J.) 

It is then natural to ask whether one could construct 
a gCFL-|-unpaired mixed phase whose components are 
electrically charged, but color neutral. This avoids the 
large color-Coulomb energy cost of mixed phases with 
colored components. Such mixed phases have recently 
been constructed in two-flavor quark matter, with un- 
paired and 2SC components [33 ■ There is still some 
color electric field (and color-charged boundary layers) 
at the interfaces where the color chemical potentials 
change rapidly as one travels from one component to 
the other, analogous to the charged boundary layers and 
ordinary electric field at the CFL /nuclear interface con- 
structed in Ref. ^, but Ref. JS^I finds that the 2SC 
-|- unpaired mixed phase does occur in two-flavor quark 
matter. However, for color neutral unpaired and gCFL 
phases, we have found that the situation is typically that 
of Fig.|5Jb): the free energy of color neutral, but electri- 
cally charged, unpaired quark matter is typically higher 
than the free energy of color- neutral gCFL, at the value 
of /ie where color-neutral gCFL is electrically neutral. 
Hence there is no value of fj,e at which oppositely charged 
phases can coexist. We have found that this is true for all 
values of Mf/fi except for a range of a few MeV just be- 
low M^/n — 130 MeV, where the neutral gCFL and neu- 
tral unpaired free energies cross. There, a mixed phase 
may arise, although as we shall discuss below it may be 
superseded by other more favorable possibilities. 

We have not eliminated all possible mixed phase con- 
structions, involving mixtures of all possible phases. 
However, over most of the gCFL regime there can be no 
mixed phase constructed from gCFL and unpaired quark 
matter. 



E. (Gapless) 2SC and 2SCus 

In this subsection, we discuss the properties of phases 
in which only two of the three flavors pair. These cannot 
compete with the CFL and gCFL phases at low values of 
M'^ / ^, but could conceivably become important at larger 
values (lower densities). 

The Fermi momenta in cold unpaired quark matter are 
ordered ppd > PFu > PFs, since the strange quark mass 
tends to decrease the strange quark Fermi momentum, 
and the down quark Fermi momentum then increases to 
preserve neutrality. Thus, the likely two-flavor pairings in 
cold three-flavor quark matter are u-d pairing (i.e. 2SC, 
with gap parameter A3 > and Ai = A2 = 0) and u- 
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FIG. 10: Gap parameters in 3-flavor quark matter for the 2SC 
phase (A3) and the 2SCus phase (A2). In each case the other 
two gaps are zero. The 2SC phase becomes gapless (g2SC) at 
M^/fj. > 113 MeV and ceases to exist at M^/fi « 130 MeV, 
and its free energy is always lower than that of unpaired quark 
matter (Fig.|^. The 2SCus phase becomes unfavored relative 
to unpaired quark matter at Mf //i > 99 MeV, and ceases to 
exist at M^ //i ~ 103 MeV, without ever becoming gapless. 



s pairing (i.e. 2SCus, with gap parameter A2 > and 
Ai = A3 = 0). 



1. Calculation 

In order to find a two-flavor pairing solution, we need 
only solve four equations (one gap and the three neu- 
trality equations). The other two gap equations are au- 
tomatically satisfied upon setting the relevant gaps to 
zero. Using the same coupling strength as in our inves- 
tigation of the gCFL phase (Aq = 25 MeV) and working 
at the same value of /i = 500 MeV, the nonzero gaps 
at M^/fi = are A3 = 31 MeV in the 2SC phase and 
A2 = 31 MeV in the 2SCus phase. As we increase M^//x, 
as long as we do not enter a gapless phase the gaps de- 
crease slowly and the simplified analysis of the 2SC and 
2SCus phases in Ref. Q should be a good guide. We 
do indeed find that our results are well approximated by 
M3 — Ms = and /ie = M^/2/i in the 2SC phase and 
fie — fi3 — /ig = in the 2SCus phase, with a free energy 
given by 
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(36) 



with Ai given by A3 in the 2SC phase and A2 in the 
2SCus phase, as predicted in Ref. pj. The free energy 
of the 2SCus phase is higher than that of the 2SC phase 
because A2 decreases more rapidly with M^/fi in the 
2SCus phase than A3 does in the 2SC phase. This cannot 
be discovered by the methods of Ref. [3] , in which these 
two phases were treated as degenerate. 
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Our results for the gap parameters are shown in Fig.llUI 
and for the free energies in Fig.|31 As in our other figures, 
we vary Ms keeping ji fixed at 500 MeV. 



2. 2SC/g2SC Results 

We find a neutral 2SC solution at low M"^ / ji^ with 
four gapped quasiparticles. At M^/fi « 113 MeV two 
of these quasiparticles become gapless, with blocking re- 
gions within which there are unpaired rd and gd quarks, 
and there is a continuous transition to the gapless 2SC 
(g2SC) phase. (Gapless 2SC was introduced in two- 
flavor quark matter in Refs. [23, 123 •) The gap pa- 
rameter then decreases rapidly until it reaches zero at 
Mg /n ~ 130 MeV and the solution ceases to exist. 

As is clear from Fig. |2| 2SC/g2SC always has lower 
free energy than unpaired quark matter, and usually has 
higher free energy than CFL/gCFL. However, we find a 
tiny window of M^/jJ. less than 1 MeV wide, very close 
to 130 MeV, in which the gapless 2SC phase has lower 
free energy than gCFL. In this regime, the one nonzero 
gap in the g2SC phase is almost zero whereas all three 
gaps are nonzero in the gCFL phase. This indicates that 
the fact that the gCFL free energy crosses that of un- 
paired quark matter almost at the same point where the 
g2SC and unpaired free energies come together is a non- 
generic feature of our model. Taken literally, our cal- 
culation predicts that as M"^ / ^ increases, gCFL is sup- 
planted by g2SC which is then almost immediately sup- 
planted by unpaired quark matter. However, treating the 
effects of Ms more accurately than we have may shut the 
tiny g2SC window completely i23|. In contrast, treating 
Ms as a chemical potential shift, as we have, but using 
Aq = 100 MeV appears to open a wide g2SC window Q, 
but this occurs in a regime where Ms ~ /i and so this 
result is not trustworthy. Also, as we discuss in the next 
section, a more general ansatz is required once one is 
at a sufficiently large Ms/fJ. that the free energy of the 
gCFL phase is close to that of unpaired quark matter, 
since there are other possible pairing patterns that likely 
become favorable. 

Finally, it is interesting to note that the 2SC solution 
in three-flavor quark matter differs from its two-flavor 
version, which requires a large /ig for neutrality given 
that there are no strange quarks present to carry nega- 
tive charge. In three-flavor quark matter, we find that at 
Ms — there is a small positive /ig and a small negative 
/ig = —fie in the 2SC phase. This happens because the 
pairing of ru, rd, gu and gd quarks increases their num- 
ber density. This contributes a positive electric charge 
and excess redness/greenness, which is compensated by 
a small positive /^e and a negative /ig. As we increase 
Ms/fi, the small /^g remains approximately unaffected 
whereas the small fj,e due to pairing is rapidly swamped 
by the larger contribution of order M^ /2// that compen- 
sates for the lack of strange quarks. 



3. 2SCus results 

We find a neutral 2SCus solution, with a gap A2 that 
decreases with increasing M^//i as shown in Fig.^| This 
solution only exists for Ms/fJ. < 103 MeV, and has a 
higher free energy than that of neutral unpaired quark 
matter for M^/fj, > 99 MeV (see Fig. El . It is always 
unfavored relative to the CFL/gCFL phase. 

It is striking that two-flavor u-s pairing, unlike the two- 
flavor u-d pairing discussed above, has no gapless phase. 
In our calculations, we find that at Ms/n > 103 MeV, 
when the 2SCus phase becomes unstable to unpairing 
(i.e. when (5/icff in the u-s sector, line 2 of table Q] be- 
comes as large as A2), there is no neutral solution with a 
smaller value of A2 (other than unpaired quark matter) . 
We do find such a "gapless 2SCus" solution in a range 
of M^/n below 103 MeV, with A2 smaller than that in 
the 2SC solution at the same Mf/fi, but the g2SCus so- 
lution is unstable: it is a local maximum of the neutral 
free energy as a function of A2. (In a figure like Fig. [7| 
this g2SCus solution would be at a local maximum of the 
solid curve.) At M^/fi = 103 MeV, the g2SCus solution 
(local maximum) meets the 2SCus solution (local mini- 
mum) at an infiection point of the neutral free energy, 
and for M^/fi > 103 MeV neither 2SCus nor g2SCus so- 
lutions exist. Unlike thegapless 2SC phase |23,|23| and 
the gapless CFL phase p, which are rendered stable by 
the constraint of neutrality, the gapless 2SCus phase re- 
mains unstable. This is presumably because u-s paired 
phases are very close to being neutral anyway (only very 
small values of fig, fi^, /ig are required to achieve neutral- 
ity _8]), so the constraint does not change the physics 
much. This can be summarized by saying that the 2SCus 
phase behaves analogously to that studied by Sarma [32| , 
even after neutrality constraints are imposed. 



IV. CONCLUDING REMARKS AND OPEN 
QUESTIONS 

The gapless CFL phase seems sufficiently well- 
motivated as a possible component of compact stars to 
warrant further study of its low energy properties and 
its phenomenological consequences: it is the phase that 
supplants the asymptotic CFL phase as a function of de- 
creasing density, and compact stars are certainly far from 
asymptotically dense. 

The low energy effective theory of the gCFL phase 
must incorporate the gapless fermionic quasiparticles 
with quadratic dispersion relations, which have number 
densities ~ fi^y/'K^T and dominate the low temperature 
specific heat, the gapless quarks with linear dispersion 
relations, with number densities ~ M^^i and the elec- 
tron excitations, with number density ~ Me-^- In con- 
trast, the (pseudo-)Goldstone bosons present in both the 
CFL and gCFL phases have number densities at most 
~ T^ . This means the gCFL phase will have very dif- 
ferent phenomenology from CFL. It will be particularly 
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interesting to compute the cooling of a compact star with 
a gCFL core, because neutrino emission wih require con- 
version between quasiparticles with hncar and quadratic 
dispersion relations. And, we expect that in a star with 
CFL, gCFL and nuclear volume fractions, the gCFL shell 
will dominate the total heat capacity and the total neu- 
trino emissivity, and thus control the (rapid) cooling. It 
will also be interesting to work out the magnetic field 
response of the gCFL phase, since the gauge boson prop- 
agators will be affected both by the gapless quasiparticles 
(all nine gauge bosons) and by the condensate (Meissner 
effect for eight out of nine.) Finally, we have left the 
study of possible meson condensation in the gapless CFL 
phase to future work. 

Although we have studied the gCFL phase in a model, 
all of the qualitative properties of this phase that we 
have focussed on appear robust. We have also offered 
a model-independent argument for the instability that 
causes the transition, and for the location of the transi- 
tion. We have used our model to show that the gCFL 
phase is favored over the two-flavor-pairing phases (2SC, 
g2SC, and 2SCus) throughout almost all of the regime 
where the gCFL phase is favored over unpaired quark 
matter. It remains a possibility, however, that the CFL 
gap is large enough that baryonic matter supplants the 
CFL phase before Mf //x > 2A. Assuming that the gCFL 
phase does replace the CFL phase, it is also possible that 
gaps are small enough that a third phase of quark matter 
could supplant the gCFL phase at still lower density, be- 
fore the transition to baryonic matter. We do not trust 
our analysis to determine this third phase. Perhaps it is 
a mixed phase of some sort, although we have ruled out 
the straightforward possibilities. Perhaps it is the gap- 
less 2SC phase |2l|,|23, as the literal application of our 
model would suggest. We should stress, in addition, that 
our model relies upon a pairing ansatz designed to study 
the instability of the CFL phase, and hence well-suited to 
the study of the gCFL phase. Determining what phase 
comes after gCFL almost certainly requires a more gen- 
eral ansatz. For example, perhaps weak pairing between 
quarks with the same flavor plays a role once gCFL is 
superseded |33, or perhaps it is the crystalline color su- 
perconducting phase ^^ |2^ jSfl 113, 113 that takes over 
from gapless CFL at lower densities. (Other possibilities 
have also been suggested |4^.') 



Recent developments |41| make the crystalline color su- 
perconducting phase look like the most viable contender 
for the "third-from-densest phase" . Previous work |4Q| 
had suggested that the face-centered-cubic crystal struc- 
ture was sufficiently favorable that its free energy could 
be competitive with that of BCS pairing over a wide 
range of parameter space, but because these indications 
came from a Ginzburg-Landau calculation pushed be- 
yond its regime of validity, quantitative results were not 
possible. The results of Ref. [4l| suggest that a crys- 
talline phase involving pairing of only two flavors is fa- 
vored over the unpaired phase by ~ 0.2^^A2g(-;/7r^ at 
Mg /jj, « 4A2SC- Here, A2SC is the gap parameter in the 
2SC phase at Ms = 0, which is 31 MeV with the param- 
eter values we have used in all our figures. This suggests 
that if we were to generalize our pairing ansatz to al- 
low the crystalline phase as a possibility, it would take 
over from gCFL at M^/fi ^120 MeV, or even some- 
what lower if the three-flavor crystalline phase, which 
no one has yet constructed, is more favorable than the 
two-flavor version. Furthermore, the authors of Ref. |4l| 
find that the crystalline phase persists until a first or- 
der crystalline— > unpaired transition at M^//i sa 7.5A2SCj 
hence over a very wide range of densities. If analysis of 
three-flavor crystalline color superconductivity supports 
these estimates, we will not have to worry about the reso- 
lution of the puzzles and possible mixed phases associated 
with the confluence of the free energies for the gCFL, 
g2SC and unpaired phases near M^/n ^ 130 MeV in 
Fig. 13 by that density the crystalline phase will already 
be robustly ensconced on the phase diagram. 
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